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Introduction 

As the connecting elements in Multi- Body structures, joints play a pivotal role 
in the overall dynamic response of these systems. Obviously, the linear stiffness 
of the joint strongly influences the system frequencies, but the joints are also 
likely to be the dominant sources of damping and nonlinearities, especially in 
aircraft and space structures. The general characteristics of such joints will be 
discussed. Then the state of the art in nonlinear joint characterization 
techniques will be surveyed. Finally, the impact that joints have on the overall 
response of structures will be evaluated. 
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Although somewhat difficult to assess, the rough order of magnitude 
of various dissipative mechanisms is shown (based on critical damping 
equaling unity). In Earth-based structures, transmission losses probably 
dominate. But in aeronautical structures, dissipation in joints begins to 
become more important. In space, in the absence of transmission losses, 
joints dominate the passive dissipation mechanisms. 


Order of Magnitude of Structural Dissipative Mechanisms 


Support transmission 
Aeroacoustic transmission 

Material damping 
Joint damping 

Active control 
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Aero 
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0 
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The potential nonlinear characteristics of a space structure are 
compared with the stiffness (normalized to unity) . In the absence of 
yielding, material nonlinearities will be on the order of fractions 
of a percent. Geometric large deflection, at least in the flexible 
modes, is small. Therefore the strong nonlinearities of the joints are 
again likely to dominate. 


Order of Magnitude of Space Structural Nonlinearities 

Material stiffness 10~^ 

Material damping 10"^ 

_2 

Geometric large deflection 10 L 

Geometric joint nonlinearity 10~^ 

Therefore joints are the largest source of passive 


DAMPING AND NONLINEARITY. 


The overall characteristics of material damping, listed below, 
coupled with the fact that the material damping is likely to be 
one-half to one order of magnitude less than joint damping, tend to 
make this a relatively less critical area in modeling. 


Material Damping Characteristics 

• Distributed with stiffness, therefore modal 

DAMPING IS PROPORTIONAL, MODES ARE REAL AND 
UNCOUPLED 

• Only weakly nonlinear, therefore approximate, 

MODELS ARE SUFFICIENTLY ACCURATE 

• Has ORIGINS IN reasonably well- understood 

MECHANISMS, E.G., THERMAL TRANSPORT, PLASTICITY 

• Is DEPENDENT ON GLOBAL FREQUENCY, AMPLITUDE 
TEMPERATURE AND HUMIDITY ENVIRONMENT 


The corresponding characteristics of joints, their discrete 
locations, strongly nonlinear behavior, and somewhat obscure micro- 
mechanics, make this a more challenging area for modeling. Despite 
frequent attempts in the history of aerospace technological develop- 
ment, no unified analysis approach to this modeling has been developed. 


Joint Characteristics 

• Not distributed, but occur at discrete locations, 

THEREFORE MODAL DAMPING IS NOT PROPORTIONAL, AND 
MODES ARE LINEARLY COUPLED AND COMPLEX 

• Strongly nonlinear, therefore modes stiffen and 

COUPLE NONLI NEARLY 

• HAS ORIGINS IN RELATIVELY POORLY UNDERSTOOD MECHANISMS, 
E . G . , MICROSLIP FRICTION, IMPACTING 


To gain some insight into this difficulty, it is useful to look 
at several proposed joint geometries for deployable space structures. 
Note that the geometries are all quite different, but all have several 
characteristics in common. There must be some amount of play in the 
joint to allow for assembly but some stiffening or locking mechanism 
to make the joint fixed when deployed. This combination of play and 
fixity leads to the impacting and nonlinear stiffness typical of such 
joints. 


Typical Joint Designs for Deployed Space Structures 




a. LaRC SNAP-JOINT UNION 


b. Rl BALL/SOCKET CONNECTOR 
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c. MIT CLUSTER SLIP-JOINT d. VOUGHT QUICK-CONNECT COUPLER 


548 







Not only do the properties of joints depend on the overall 
geometry of their design, but these properties depend on a number 
of details. The surface of the contacting elements can depend, 
for example, on the quality of machining, the load and wear history, 
and the duration on orbit. Even nominally identical joints can have 
a statistical variation due to manufacturing tolerance. Therefore, 
in realistic assemblies, direct calculation of properties is somewhat 
unproductive. 


• Joint properties depend on very local details 

1. SURFACE FINISH, LUBRICATION, OUTGASSING 
AND OXIDATION 

2. WEAR AND TRIBOLOGY 

3. PRECISION OF FIT AND ALIGNMENT 
q . PRELOAD AND INITIAL DEFORMATION 
5. LOCAL THERMAL DEFORMATIONS 

• Joints of identical materials can have very 

DIFFERENT BEHAVIOR 

§ Nominally identical joints may have a statistical 

VARIATION IN BEHAVIOR 

Therefore the detailed calculation of joint characteristics 

FROM FIRST PRINCIPLES IS UNPRODUCTIVE, 
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A more common approach to characterization is a hybrid of simplified 
modeling and experimentation, A set of experiments is run, yielding some 
data on the force transmission of the joint. Concurrently, several 
postulated models of the joint are developed. Often this is somewhat 
interactive, i.e., after the data are evaluated, refined models are 
postulated. The force characteristics, or structural response of the 
postulated model, is then compared with the experimental data, and 
some fit of the model to data is performed. Based on this fit, the 
parameters of the model are available for use in the overall structural 
model. 


ALTERNATIVE TO DIRECT CALCULATION 
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A number of different models of joint behavior can be postulated. 

Three of the common ones are shown below. The first is Coulomb friction, 
in which the normal force, and therefore the frictional drag force, remains 
constant. In displacement dependent friction, the normal and frictional 
drag forces vary with displacement. This model is probably more realistic 
for jointed trusses in the absence of thermal and gravity loads. When the 
deadband closes, impacting occurs, and a sharp jump in damping and stiff- 
ness (not shown) occurs. 


Example Postulated Piecewise Linear Models 
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The principal characteristics of four procedures for identification 
of nonlinear elements are shown. The first two are extensions of tech- 
niques developed for linear systems and are more easily extendable to 
multi-dof-models . However, they are probably only appropriate for weak 
nonlinearities. The latter two are currently limited to single-dof systems, 
but can handle stronger nonlinearities. A more detailed explanation of 
each will follow. 


Experimental Identification Procedures 


Procedure 

Measurement 

Frequency Domain 

Modal Identification [2] 
(Ewins,Imp,Col.) 

F, x vs “ 

Transient Time 

Domain t3i 
(Horta, Juang, LaRC) 

F, X VS T 

Classical Force- 

Displacement mi 
(Son i, Udri) 

F vs x 

Force-State (5] 
(Crawley, MIT) 

F vs x, x 
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The extended frequency domain modal identification procedure was 
developed simply to uncover the presence of nonlinearities in modal 
test data. Therefore, the output is limited to indications of the 
presence, strength and type of nonlinearity. It is best used as a 
diagnostic tool in checking the consistency of test data. 


Extended Frequency Domain Modal Procedure t2i 


• Postulate nonlinearity and calculate loss factor 

USING NyQUIST PLANE RESPONSE 

• Measure response at resonance and calculate loss 

FACTOR 

• If loss factor is inconsistent (i.e., not constant), 

CHOOSE POSTULATED NONLINEARITY WHICH BEST FITS 
OBSERVED BEHAVIOR 

• Output - approximate indication of type and degree 

OF NONLINEARITY IN MODAL RESONANCE. 


The technique is a direct extension of the procedure for extracting 
frequency and loss factor parameters from transfer functions, as presented 
in the Nyquist plane (b). A simple, single-dof response appears as a 
perfect circle in this representation, tangent to the real axis at the 
origin. Any deviation from this circle is due to a nonlinearity (or pre- 
sence of multiple poles). The loss factor (damping ratio) can be calculated 
by choosing pairs of points about U) , forming a matrix of computed values. 
Figure (c) is a graphical representation of loss factor calculated on this 
matrix. For a linear system, this surface would be flat. The shape shown 
is typical of a system with Coulomb friction. 


CALCULATED FREQUENCY RESPONSE (A) AND NYQUIST REPRESENTATION (B) FOR A SDOF 
SYSTEM WITH COULOMB FRICTION. THE LOSS FACTOR (C) IS INFERRED FOR A RANGE 
OF FREQUENCY SPREAD ABOUT THE RESONANCE. 
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Likewise, the existing time domain techniques are extensions of 
techniques developed for linear systems. These techniques generally 
examine the transient free response to extract system mode shapes and 
frequencies. Weak nonlinearities appear as a frequency with a number 
of higher harmonics. Each type of nonlinearity has such a signature. 


Extended Eigensystem Realization Algorithm t3i 

• Postulate nonlinearity and calculate Fourier 

CONTENT OF TRANSIENT FREE RESPONSE 

• Measure free response and identify Fourier 
CONTENT WITH ERA 

• Compare measured higher harmonic content of 

of modal response with signatures of postulated 
nonlinearities 

• Output - approximate indication of type and degree 
of nonlinearity in response, 



Four example cases are shown, all typical of a stiffening or softening 
spring. The Fourier components of the free response of such a spring in a 
spring mass system were calculated. The calculated response was also fed 
as simulated data to the ERA program and the harmonics of the response 
calculated. Good capability to reconstruct the signature of a known non- 
linearity is shown. However, the recognition of the signature of an unknown 
nonlinearity is still under development. 


Four generic nonlinear joints (a) and 
the Fourier content of their transient 

DECAY FROM ANALYSIS AND ERA IDENTIFICATION 
OF COMPUTED RESPONSE (b) . 



(A) 



FREQUENCY Hz . 

COMPONENT AMP. 

CASE NO. 

ERA 

ANALYSIS 

ERA 

ANALYSIS 


0. 135 

0. 134 

0.986 

0. 987 

1 

0.404 

0. 403 

0.015 

0.015 


0.673 

0.672 

-0. 002 

-0.002 


0. 125 

0. 125 

0.979 

0.979 

o 

0.375 

0.375 

0.023 

0. 023 


0.625 

0.625 

-0.003 

-0.003 


0.096 

0.096 

1 . 030 

1 . 030 

3 

0. 289 

0.289 

-0.035 

-0.035 


0.482 

0. 482 

0. 006 

0. 006 


0.456 

0.456 

0.961 

0.961 

4 

1 . 369 

1 . 369 

0.038 

0.038 


2.281 

2.281 

0.001 

0.001 


(b) 
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The classical approach to the problem is, of course, to simply put a 
joint in a testing machine and develop force-stroke data, as shown on the 
next page. From such data secant modulus and average loss factor can be 
calculated. The limitation is that such properties are already smeared, 
or averaged over the stroke, and no tangent or point properties can be 
determined. Further, the dependence of the force on the rate of change 
of stroke is lost. 


Quasi-Steady Force Deflection Procedure UJ 

• Postulate nonlinearity and calculate its F vs x 

BEHAVIOR 

• Measure F vs x behavior for one x max and and 

CALCULATE EFFECTIVE STIFFNESS AND LOSS FACTOR 

• Repeat at different amplitudes and frequencies 

• OUTPUT - EFFECTIVE STIFFNESS AND LOSS AS A 
FUNCTION OF FREQUENCY AND AMPLITUDE. 
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LOAD 


Typical Force-Stroke cycles of the 

THREE STRESS RESULTANTS OF A JOINT 
AT ONE LOAD AND FREQUENCY, 



AXIAL TENSION 


TRANSVERSE BENDING 


INPLANE BENDING 



The force-state mapping procedure is designed specifically to identify 
strong nonlinearities in joints and addresses the two limitations of the 
classic Force-Stroke measurement. The dependence of transmitted force on 
both displacement and velocity is explicitly determined, and local or tangent 
values produced* At the current time the procedure is limited to single-dof 
systems . 


Force-State Happing Procedure 

• Postulate nonlinearity and calculate 
F vs x, x behavior 

/ 

• Measure F vs x, x over expected range 

• Fit postulated surfaces to date in 
Force-State space 

• Output: 

1. Raw data for data look-up 

2. Global fit parameters for analytic 

REPRESENTATION 

3. Local equivalent secant moduli 

FOR LINEARIZED REPRESENTATION 
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The force-state maps of three simple systems are shown. A spring (a) 
produces a plane with a slope against x, but no change in x. A linear 
viscous damper would produce a plane with a slope against x, and no 
change in x. Thus any linear element (i.e., spring and damper) will have 
a map which is a flat plane in force-state space. Any deviation from a 
plane is indicative of a nonlinearity. 

Two common nonlinearities are shown in (b) and (c). The cubic spring 
nature is clear in fig. (b). Figure (c) shows the map of Coulomb friction, 
which is independent of x, and takes on the sign of the velocity. 


Force-State Maps of: 

a) linear spring 

b) cubic spring 

c) Coulomb friction 



i 
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(c) 


ORIGINAL PAGE IS 
OF POOR QUALITY 



The force-state maps of a real joint are shown below. The joint is 
of a quick-disconnect-pin and clevis type, similar to the Vought connector 
shown earlier. The figures on the left (a and b) show the characteristic 
without a stiffening sleeve. Note the step at x equals zero, indicative 
of friction. With the addition of a sleeve, the joint becomes stiffer 
(note the change in vertical scale), and the cubic stiffness of the sleeve 
begins to dominate. A strong dissipative nature is still obvious from the 
hysteresis loops in fig. (d). 


The Force-State and Force-Displacement 

CURVES OF A CLEVIS-PIN-TYPE CONNECTOR 
WITHOUT (A AND B) AND WITH (C AND D) 

A REINFORCING SLEEVE. 


30. 6 


250.4 



In an effort to fit a postulated model to the data on the previous 
page (i.e., the joint with sleeve), three successively refined models 
were used. In fig. (b) , a cubic, plus linear, spring term approximately 
matches the shape but, of course, has no dissipative nature. The addition 
of friction introduces the classic hysterectic step. Finally, with the 
introduction of linear damping, the measured data are closely reproduced. 


Successive approximations of actual data (a) by a cubic spring (b), 

CUBIC SPRING, PLUS FRICTION (c) AND CUBIC SPRING, FRICTION AND LINEAR 
DAMPING (d). 






WSPUCniENT (U) 







The requirements for efficient computation place several requirements 
on the identification scheme. It is highly desirable to have available 
the force-state information for direct pseudo-force computation. 


Computational Considerations 1 63 

• Three computational approaches to including 

THE JOINT NONLINEARITY CAN BE CONSIDERED 

1. Homogeneous nonlinearity, explicit operator 

2. Homogeneous nonlinearity, implicit operator 

3. Pseudo-force representation [ 7 ] 

§ In all three, but especially in the pseudo-force 

METHOD, IT IS NECESSARY TO HAVE THE JOINT CHARAC- 
TERISTIC IN TERMS OF JOINT STATE VARIABLES. 

• If ONLY AVERAGE, OR SECANT PROPERTIES ARE KNOWN, 
THEN CONSIDERABLE ITERATION IS REQUIRED, AND 
TRANSIENT ANALYSIS MAY NOT BE ACCURATE. 
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As an example problem, a four-bay truss, connected by joints, is 
modeled as a four segment beam, pinned in translation. In rotation it 
is constrained by a linear spring and damper. 


Effects of Joints on Modal Properties 
Model a connected 2-D truss 



AS A PINNED BEAM OF 4 ELEMENTS WITH 
ROTARY SPRINGS AND DAMPERS. 
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When the natural frequencies of the system are plotted versus 
nondimensional joint stiffness, their trends are apparent. Of course, 
all modes stiffen as k is increased- Some modes , such as #4, are only 
slightly affected, while others, such as #7, are strongly affected. 

The lowest eight modes are asymptotic to a constant frequency, while the 
highest three continue to rise as k increases. 

Stiffening effects of joints as a function 
0F k joint* 


FREQUENCY vs JOINT STIFFNESS 

(4— Element Beom with Pin-Joints) 




The addition of linear joint damping has some surprising results. 
Note that in only three modes, 7, 10 and 11, is the damping roughly 
proportional, i.e., the pole is driven to the real axis. In most modes, 
the root damps, then asymptotically stiffens and loses damping. In one 
mode, //9, the frequency drops. 


Locus OF ROOTS FOR INCREASING LINEAR 
JOINT DAMPING, FOR K lrtIkl -r = 0.3 EIA 





Finally, this figure shows an interesting application of the force- 
state map to Earth testing of space structures. Suppose a structure was 
suspended in one gravity in such a way that the gravity load caused a 
steady deflection. The small displacement vibration would then take place 
about this "Earth IC," and would have the effective stiffness and damping 
shown. In space, in the absence of gravity loads, there would be no steady 
deflection and the effective K and C would be about a "Space IC," as shown. 
For a generally nonlinear joint, these properties could be completely dif- 
ferent from those of the Earth test, leading to differences in dynamic 
behavior on orbit when compared to those measured on Earth. 


Use of the force-state map to determine 

THE EFFECTIVE STIFFNESS AND DAMPING IN A 
JOINTED STRUCTURE, AS WOULD BE MEASURED 
ON EARTH AND IN SPACE. 


F ~MA EFFECTIVE K AND C 
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Summary 


• Detailed modeling of micromechanics 

OF THE JOINT NOT PRODUCTIVE 

• Development of simple generic models 

USEFUL 

• Improved nonlinear identification 


NECESSARY. 
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